Dielectric properties of a two dimensional binary sytem with ellipse inclusions 
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A two-dimensional binary composite system composed of ellipse inclusions and a host medium 
is considered. Dielectric permittivity of the system is calculated as a function of orientation angle, 
the volume fraction of the inclusions and their excentricity using the finite element method. It was 
observed that both the orientation of the inclusions in the field and their excentricity have significant 
effects on the dielectric permittivity. 



Binary systems composed of two media have been stud- 
ied through the history due to their simplicity. Dielectric 
properties of such systems show different behaivour de- 
ending on the distribution and shape of the inclusions 
Jj. Materials used in electrical insulation sytems are 
usually composed of several materials. Therefore, behav- 
ior of such systems in electrical fields are crucial for the 
insulation. 

In this paper, we determine the effective parameters 
of a two diemsional structure composed of a main phase 
and ellipse inclusions. We have assumed that the struc- 
ture is periodic where the inclusions sit on lattice points 
of a square array. With such an assumption the sys- 
tem can be subdivided into smaller units i.e. square 
cells. This approach was first introduced by Rayleigh 0] 
and it has been used in numerical and analytical tech- 
niques j|, 0. In this model, inclusions were arranged 
on rectangular lattice points. One can assume cylinder- 
ical inclusions in two-dimensions and pherical inclusions 
in three dimesions. Analytical solutions for two compo- 
nent and three component systems with different inclu- 
sion geometries (infintely long cylinders) were given in 
Refs. H 0, || . Moreover, frequency dependent charac- 
teristics, lattice effects and random structures were pre- 
sented in Refs. ]|, g|. 

Mean field theory approaches and solutions for fre- 
quency dependent effective medium properties of a mix- 
ture with ellipsoidal inclusions were studied by Sillars 
PJ and Steeman and Maurer pi. The former one as- 
sumed that the inclusions were much more conductive 
than the matrix material, and the later one introduced 
the conductivity of the matrix. In both cases they have 
assumed that the ellipsoidals were parallel or perpendic- 
ular to the applied electric field. Therefore, we concen- 
trated on the orientation of the ellipses in the periodic 
structure. First, all the ellipses were assumed to have 
the same orientation. Later, we change the orientation 
of the neighbor ellipse which made the system doubly pe- 
riodic. We exclude the effects of ellipse concentration on 
effective parameters in the mixture. 

Let ellipses with dielectric permittivities Ei be located 
in an unbounded medium with e m . The ellipses are 
placed at the square lattice sites as presented in Fig. 
pi The solution of the Maxwell equations without any 
magnetic field presence, will assist to estimate the ef- 
fective properties of the mixture. Assumptions of the 
ellipse in two dimensions and the ellipsoidal in three di- 
mensions can explain effective properties of composite 
structures better than considerations with spherical or 



FIG. 1: Unit cell. 



cylinderical inclusions, since real composites have arbi- 
trary shaped inclusions. Moreover, the interaction be- 
tween ellipse disks will be different depending on the ori- 
entation of the neighboring ellipses. 

Unit ellipse in xy-coordinates is given by (Fig. 0) 



(x-xp (y-y ) 2 
a 2 b 2 



(1) 



where xq and j/n are the center coordinates of the ellipse, 
and a and b are the maximum radii in x and y directions. 
One can also use polar-coordinates, 



(2) 



Moreover, if one includes the orientation of the ellipse, 
angle from the x-axis /3, in the matrix medium, the x and 
y points will be multiplied by the transformation matrix, 
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A point, L(x, y) can be calculated in this way. The shape 
factor, eccentricity (e), of the ellipse is 
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FIG. 2: Geometrical parameters for ellipses. 
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FIG. 3: Equi-concentration lines. 



intrinsic electrical properties i.e. dc conductivity, a and 
polarization, e of the medium. Excluding the dipolar re- 
laxations in the material, the complex dielectric constant 
becomes 



EqUJ 



(6) 



When two materials form an interface that satisfy the 
condition e\/a\ =^ £2/02, a relaxation process due to 
the interfacial (Maxwell- Wagner-Sillars) polarization is 
observed. Since we have excluded the other relaxation 
processes coming from the materials, the dielectric per- 
mittivity of the mixture will be swithed after some fre- 
quency. It will be the same for the conductivity. 

The calculation of the effective properties of the mix- 
ture from the known material properties of the con- 
stituents are electromagnetic problem which involve the 
solution of Maxwell equations with boundary conditions. 
Numerical solutions of electrostatic problems within a 
non-conducting medium are based on the solutions of 
Poisson's equation 



V • (eV$) = -Q 



(7) 



where e, $ and Q are the electrical potential, permit- 
tivity of the medium and total charge in the considered 
region, respectively. However, if the medium has conduc- 
tive regions, where no free charges and sources of charges 
are allowed, then, the solutions are given by 



V • (<rV$) = 



(8) 



where a is the conductivity of the region. When the 
medium is a mixture of these two cases, it consists of 
capacitive and resistive components. The solution is then 
given by a time dependent and complex electric potential 
in the region with the coupling of Eqs. (0) and (g). 



V • (crV$) + V • 
V-((<7 
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(9) 



In this equation, the boundaries of the unit cell as well 
as the area fraction of the inclusion ellipse are the con- 
straints for the values of a and b. Since, the ellipse in- 
clusions are hard disk-like structures, they are not al- 
lowed to overlap. This leads to a < 1/2 and b < 1.2 in 
the unit cell. The resulting a and b values versus equi- 
concentration lines and the eccentricity are presented in 
Fig. H and Fig. [|. However, for orientation angles, (3 
other than [3 7^ ri7r/2, n being an integer, the previous 
constraints will be different, which is not included in this 
paper. The area fraction, q, of ellipse in the unit square 
medium is 



q = nab 



(5) 



To study the frequency dependent effective properties, 
for example, effective dielectric permittivity and conduc- 
tivity, of the mixture, finite element method (FEM) nu- 
merical technique was used in two dimensions. The com- 
plex dielectric constant, e* carries information about the 




FIG. 4: Eccentricity as a function of a and b. 
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FIG. 5: Complex dielectric susceptibilities, x*i °f a binary 
system with different ellipse orientations and eccentricities at 
q = 0.1. The angle for neighboring ellipses is the same; (a) 
/?i,a = n /2 and (b) 0\ y % — 0. The labels I and II represent x' 
and ^''-values for e = 0.980. In a similar way, the labels Ml 
and IV represent x' and x" _v alues for e = 0.965. The labels 
V and VI represent \ an d x" _vames f° r a circular inclusion, 
e = 0. 
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FIG. 6: Complex dielectric susceptibilities, x*i °f a binary 
system with different inclusion orientations at q — 0.1 when 
e = 0.980. The angle for neighboring ellipses is the same 
(a) fa = 2 = tt/2 (labels I and II) and t = 2 = (labels 
III and IV), and (b) ft = #2 = tt/3 (labels III and IV) and 
Pi=fo = 7r/6 (labels III and IV). The labels V and VI are for 
a circular inclusion, e = 0. The odd and even labels represent 
x' Mid x"- va l ues > respectively. 



where i is y/—l, no free charges are allowed in the region 
due to resistive component. Last two equations are the 
continuity and the Laplace equations for a medium with 
a complex dielectric constant. 

Finite element method was used in order to solve the 
equations with the help of a field calculation software pT[. 
In this method, the region is meshed using an adaptive 
meshing technique with triangles and then the equations 
are solved for each triangle by interpolating the potential 
and its normal derivative (electric field) using boundary 
conditions |T^| . In order to have a reliable numerical so- 
lution, linear, quadratic and cubic solution sets are com- 
pared for different mesh sizes and results of the cubic 
solutions are used when there were not any differences in 
the solutions. The boundary conditions were chosen due 
to symmetry of the unit cell, which were repeated in one 
direction cr(d&/dn) = —j, where j is the current den- 
sity, and rms voltage difference, AV, of 1 V were applied 



on the other boundaries, respectively. It is illustrated in 
Fig. 0. 

Frequency dependent capacitance was calculated using 
the total charges, Qt per unit length on one of the voltage 
applied boundaries by the Gauss' law, 



C(f) 



QtU) 

Ay 



D- n dl 



(10) 



where, D is the obtained electrical displacement normal 
to the line incriment dl of the region, I |13|] , and it is 
equal to, 



D = -e 



<9$ 
On 



(11) 



Finally, frequency dependent real part of the relative 
permittivity, e'(/), was obtained by this value over the 



geometrical capacitance of the primitive cell, Co 



e'(f) = 



C(f) 
Co 



(12) 



In the same way, frequency dependent conductance, 
G(f) was calculated using the total normal current, /„ 
per unit length on one of the voltage applied boundaries 
by Ampere's law, 



G(f) 
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(13) 



where, I is the obtained current to the line incriment d\ 
of the region, I, and it is equal to, 



<9$ 

on 



(14) 



Frequency dependent imaginary part of the dielectric 
permittivity, e"{f) is then 



e"{f) 



G(f) 
27re / 



(15) 



In the simulations, the material parameters of the me- 
dia were s± — 4 and £2 = 40, and conductivities <J\ = 
lCT 12 frW 1 and o x = lO" 10 fl' 1 m~ 1 . The boundary 
conditions were assigned such that when j3x,2 = tt/2 the 
inclusions were parallel to the electric held, E. 

The influence of shape of inclusions were investigated 
by considering two e- values at q = 0.1. Moreover, 
circular-inclusions, e = 0, were also considered for com- 
parison. In Fig. a, the calculated complex dielectric 
susceptibilities, x* , are presented for e = 0.980 and 
e = 0.965 while fix and ^-values were kept constant, 
/h,2 = t/2 and 0x,2 = 0, respectively. When the e-value 
was close to 1, the shape of the inclusion was needle- 
like. Moreover, when the orientation of the inclusions 
were parallel to E, (3x,2 = tt/2, they had not only higher 
X'- values, but also higher losses, x" \ than the circular- 
shape. If the x* -values of two elliptical-shapes were com- 
pared, the ellipse with higher e-value had the highest 
X*-values. If the opposite case was taken into consid- 
eration, in which the orientations of the inclusions were 
perpendicular to E, (3x,2 = 0, then, both shapes had sim- 
ilar x'-curves which were lower than that of the circular- 
shape. The x"-values of the all considered cases showed 
the ohmic losses which hid the losses due to interfacial 
polarization. 

In Fig. m, a more general case, when j3x 7^ P2, is pre- 
sented. When Px — 0, the influence of the orientational 
state, @2 , of the neighboring inclusions was affecting the 
X*-values. The effect was, on the contrary, not that sig- 
nificant when (3x — 7r/2, as illustrated in Fig. |7|. Finally, 
the influence of orientation of inclusions on dielectric per- 
mittivity at high frequencies, Ehf was considered, and is 
displayed in Fig. H. The £/j/-values increased when the 



orientations of nearest neighbors were close to the direc- 
tion of the applied field, E. 

In order to investigate the shape and orientation of 
inclusions, elliptical inclusions were considered. Needle- 
like inclusions had higher dielectric strengths, Ae, and 
observable dielectric losses, x" > compared to the circular 
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FIG. 7: Complex dielectric susceptibilities, x*i °f a binary 
system with different orientations at q = 0.1 when e = 0.980. 
The orientation of one ellipse is kept constant (a) /3i = and 
(b) Px — 7r/2, and the orientations of the nearest neighbors, 
,02 are altered. The angle values are /32 = [0, 7r/6, 7r/3, 7r/2] 
for the labels. The odd and even labels represent %' and x"~ 
values, respectively. 



shapes. The effects were larger when the orientation of 
the longer-axis of the inclusions were parallel to the elec- 
tric field. Moreover, when the orientations of the near- 
est neighbors were considered, higher dielectric strengths 
were observed for neighbors oriented close in the direc- 
tion of the applied field. 
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FIG. 8: Dielectric permittivity at high frequencies, euf, of 
a binary system with different orientations of inclusions at 
q = 0.1 when e = 0.980. The /3i-value is kept constant and 
/fe-values are altered. The labels I, II, III and IV represent the 
cases when ft = [0, 7r/6, 7r/3, tt/2], respectively. The solid- line 
(V) is the £h/-value for a circular-inclusion. 
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